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ABSTRACT 

We present an analysis of the chemical abundances of the star Tycho G in the 
direction of the remnant of supernova (SN) 1572, based on Keck high-resolution 
optical spectra. The stellar parameters of this star are found to be those of a 
G-type subgiant with T cS = 5900 ± 100 K, log(#/cm s2) = 3.85 ± 0.30 dex, and 
[Fe/H] = —0.05 ± 0.09. This determination agrees with the stellar parameters 
derived for the star in a previous survey for the possible companion star of SN 
1572 (Ruiz-Lapuente et al. 2004). The chemical abundances follow the Galactic 
trends, except for Ni, which is overabundant relative to Fe, [Ni/Fe] = 0.16 ± 0.04. 
Co is slightly overabundant (at a low significance level). These enhancements in 
Fe-peak elements could have originated from pollution by the supernova ejecta. 
We find a surprisingly high Li abundance for a star that has evolved away from 
the main sequence. We discuss these findings in the context of companion stars 
of super novae. 

Subject headings: stars: abundances — supernovae: general — stars: evolution 
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1. Introduction 

Type la supernovae (SNe la) are the 
best known cosmological distance indica- 
tors at high redshifts. Their use led to the 
discovery of the currently accelerating ex- 
pansion of the Universe (Riess et al. 1998; 
Perlmutter et al. 1999); see Filippenko 
(2005) for a review. They were also used 
to reveal an early era of deceleration, up 
through about 9 billion years after the big 
bang (Riess et al. 2004, 2007). Larger, 
higher- quality samples of SNe la, together 
with other data, are now providing increas- 
ingly accurate and precise measurements of 
the dark energy equation-of-state parame- 
ter, w = P/(pc2) (e.g., Astier et al. 2006; 
Riess et al. 2007; Wood-Vasey et al. 2007; 
Kowalski et al. 2008). 

Though the increase in the empirical 
knowledge of SNe la has led to an enor- 
mous advance in their cosmological use, 
the understanding of the explosion mecha- 
nism still requires careful evaluation (e.g., 
Hillebrandt & Niemeyer 2000; Wheeler 
2007). While in Type II SNe we have the 
advantage that the explosion leaves a com- 
pact star to which surviving companions 
(if they exist) often remain bound, thus 
enabling a large number of studies (e.g., 
Martin et al. 1992; Israelian et al. 1999), 
in SNe la the explosion almost certainly 
does not produce any compact object. 

To investigate how the explosion takes 
place, we examine the rates of SNe la at 
high redshifts, or we can take a more direct 
approach and survey the field of histori- 
cal SNe la. The latter strategy has been 
followed since 1997 in a collaboration that 
used the observatories at La Palma, Lick, 
and Keck (Ruiz-Lapuente et al. 2004). 



The stars appearing within the 15% in- 
nermost area of the remnant of SN 1572 
(0.65') were observed both photometrically 
and spectroscopically at multiple epochs 
over seven years. The proper motions of 
these stars were measured as well using im- 
ages obtained with the WFPC2 on board 
the Hubble Space Telescope (HST) (GO- 
9729). The results revealed that many 
properties of any surviving companion star 
of SN 1572 were unlike those predicted by 
hydrodynamical models. For example, the 
surviving companion star (if present) could 
not be an overluminous object nor a blue 
star, as there were none in that field. Red 
giants and He stars were also discarded as 
possible companions. 

A subgiant star (G2 IV) with metal- 
licity close to solar and moving at high 
speed for its distance was suggested as the 
likely surviving companion of the explod- 
ing white dwarf (WD) that produced SN 
1572 (Ruiz-Lapuente et al. 2004). This 
star, denoted Tycho G, has coordinates 
a = 00 h 25 m 23.7 s and 5 = +64 o 08'02" 
(J2000.0). Comparisons with a Galactic 
model showed that the probability of find- 
ing a rapidly moving subgiant with solar 
metallicity at a location compatible with 
the distance of SN 1572 was very low. 
A viable scenario for the Tycho SN 1572 
progenitor would be a system resembling 
the recurrent nova U Scorpii. This sys- 
tem contains a white dwarf of Mwd = 
1.55 ± 0.24 M and a secondary star with 
M 2 = 0.88 ± 0.17 M orbiting with a pe- 
riod P orb ~ 1.23 d (Thoroughgood et al. 
2001). 

On the other hand, based on analysis of 
low-resolution spectra, Ihara et al. (2007) 
recently claimed that the spectral type of 
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Tycho G is not G2 IV as found by Ruiz- 
Lapuente et al. (2004), but rather F8 V. 
In addition, Ihara et al. (2007) suggested 
the star Tycho E as the companion star of 
the Tycho SN 1572, but we consider their 
conclusion to be unjustified. They base it 
on a single Fe I feature at 3720 A in a spec- 
trum having a signal-to-noise ratio (SNR) 
of only ~ 13. The short spectral range of 
their data (3600-4200 A) and the low res- 
olution (A/AA ~ 400, with a dispersion of 
5 A/pixel) add to our concern. 

The present work resolves questions re- 
garding the metallicity and spectral clas- 
sification of Tycho G. We concentrate on 
providing both the chemical composition 
of the star and the stellar parameters as 
derived from new data. In addition, we 
present a detailed chemical analysis of Ty- 
cho G aimed at examining the possible pol- 
lution of the companion star by the su- 
pernova. Although more data are still re- 
quired to definitively settle this issue, the 
overall properties of this star remain con- 
sistent with being the surviving companion 
of SN 1572. 

2. Observations 

In order to perform an abundance anal- 
ysis of Tycho G, spectra of it were ob- 
tained with the High Resolution Echelle 
Spectrometer (HIRES; Vogt et al. 1994) 
on the Keck I 10-m telescope (Mauna Kea, 
Hawaii). Nine spectra (with individual ex- 
posure times of 1200 s and 1800 s) were 
obtained on 10 September 2006 (UT dates 
are used throughout this paper) and four 
spectra (with exposure times of 1800 s) 
on 11 October 2006, covering the spectral 
regions 3930-5330 A, 5380-6920 A, and 
6980-8560 A at resolving power A/AA ~ 



50,000. 

The spectra were reduced in a stan- 
dard manner using the makee package. 
We checked in each individual spectrum 
the accuracy of the wavelength calibra- 
tion using the [O I] A6300.3 and A5577.4 
night-sky lines and found it to be within 
0.2 km s _1 . After putting each individ- 
ual spectrum into the heliocentric frame, 
we combined all of the spectra from each 
night. Then, each night's final spectrum 
was cross-correlated with the solar spec- 
trum (Kurucz et al. 1984) properly broad- 
ened with the instrumental resolution of 
~ 6.6 kms" 1 . The final HIRES spectrum 
of Tycho G has a SNR of ~ 20 at 5200 A, 
~ 30 at 6500 A, and ~ 50 at 7800 A. This 
spectrum was not corrected from telluric 
lines since it was unnecessary for the chem- 
ical analysis (see § 13. 3p . 

With the aim of determining the spec- 
tral classification of several targets in the 
Tycho SN 1572 field, on 12 December 2007 
we obtained spectra with the Low Resolu- 
tion Imaging Spectrometer (LRIS; Oke et 
al. 1995) on the Keck I telescope. We used 
the 400/3400 grism, the 400/8500 grat- 
ing, and the D560 dichroic, covering the 
range 3180-9150 A. The full width at half- 
maximum (FWHM) resolution was 6.4 A 
in the red part (A > 5500 A) and 6.3 A 
in the blue, with respective dispersions of 
1.86 A pixel" 1 and 1.09 A pixel" 1 . Ty- 
cho E, F, G, and D were observed in single 
exposures of 700, 200, 400, and 450 s, re- 
spectively; see Ruiz-Lapuente et al. (2004) 
for star identifications. The SNR of the 
individual spectra is 21, 37, 31, and 25 
at 4600 A for Tycho E, F, G, and D, re- 
spectively. These spectra were reduced us- 
ing standard techniques (e.g., Foley et al. 
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Fig. 1. — Left panel: Keck LRIS spectrum of Tycho G compared with three template spectra 
having spectral types G2 IV, F8 I, and F6 V taken from a library of low-resolution stellar 
spectra (Jacoby et al. 1984). The absorption feature near 6900 A in the template spectra 
is the telluric B band. Another, weaker telluric feature is visible near 7200 A. Right panel: 
LRIS spectrum of Tycho G dereddened with E(B — V) = 0.66, 0.76, and 0.86 mag, compared 
with spectra from the low-resolution library of Jacoby et al. (1984). The best fits to the 
different dereddened spectra of Tycho G are also shown. 



2003), including removal of telluric absorp- 
tion lines. 

The spectrum of Tycho G is shown in 
Figure 1; spectra and classifications of the 
other stars are available in Appendix [A] 
In general, our results do not support the 
spectral classifications suggested by Ihara 
et al. (2007), which are based on low- 
resolution (A/ A A « 400), low-SNR spec- 
tra obtained with the instrument FOCAS 
at the Subaru Telescope over a short spec- 
tral range (3600-4200 A). Moreover, their 
relative flux calibrations are questionable 
since they apply a correction to the slope 
of each spectrum that is fitted when com- 
paring with a template. Thus, we consider 
their spectral classifications to be unreli- 
able. 



3. Stellar Parameters 

3.1. Low- Resolution Spectra 

The blue and red parts of the LRIS spec- 
tra were merged and then rebinned at a 
scale of 2 A pixel -1 . We compared the 
spectra of the four targets with a library of 
low-resolution stellar spectra from Jacoby, 
Hunter, & Christian (1984). We chose this 
library because the spectra have a similar 
resolution (~4.5 A), with a dispersion of 
1.4 A pixel -1 , and cover a similar spectral 
region (3510-7427 A) as our LRIS obser- 
vations. 

The low-resolution comparison with 
template spectra depends on the color 
excess applied to the LRIS spectra. We 
dereddened the calibrated spectra accord- 
ing to the parameterization of Cardelli, 
Clayton, & Mathis (1989), including the 
update for the near-ultraviolet given by 
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O'Donnell (1994). Moreover, there is un- 
certainty in the spectral classification of 
each template. Thus, the spectral classi- 
fication of stars using low-resolution spec- 
tra must always be considered only as an 
initial guess. The best fit is found for a 
reddening E(B — V) = 0.76 mag, with 
X2 U = 3.68 (see below for details). This 
value of reddening was used in the com- 
parison shown in Figure [TJ In the library 
of Jacoby et al. (1984), there are not many 
spectra with luminosity class IV (only F0, 
F3, G2, G5 IV), whereas luminosity class 
V is very well sampled with several tem- 
plates for each spectral type. If only low- 
resolution spectra are used to derive the 
spectral and luminosity classes, the results 
may be erroneous, since they depend on 
reddening and on the flux calibration of 
the observed spectra. 

Nonetheless, clS 8b first step, we com- 
pare the LRIS spectrum of Tycho G with 
the low-resolution spectra of three tem- 
plate stars having different spectral types 
(Fig. CD). To evaluate the goodness of fit, 
we employ a reduced x2 statistic, 

where iV is the number of wavelength 
points, M is the number of free param- 
eters (here one), fi is the template flux, Fi 
is the target flux (in this case, Tycho G), 
and o"j = 1/SNR. The SNR was estimated 
as a constant average value in the contin- 
uum over the range 4610-4630 A; this re- 
gion was also used to normalize both target 
and template spectra. The spectral regions 
used to evaluate %2 are 3850-6800 A and 
7000-7400 A, avoiding strong telluric lines 
near 6900 A present only in the template 



spectra. 

The best fit is found for the G2 IV tem- 
plate with x^v = 3.68, whereas the F8 I 
and F6 V templates provide x^u = 6.39 
and 9.66, respectively. Ihara et al. (2007) 
have claimed that the spectral type of Ty- 
cho G could be either F8 I, F8 V, or 
F6 V, rather than G2 IV as found by Ruiz- 
Lapuente et al. (2004). However, our anal- 
ysis shows that those other spectral types 
provide worse fits than the G2 IV template. 

In the following sections we use more ac- 
curate classification methods, which are in- 
dependent of reddening, to derive not only 
the spectral type of Tycho G but actually 
the effective temperature and surface grav- 
ity, using a high-resolution spectrum. 

3.2. Ha Profile 

The wings of Ha are a very good tem- 
perature indicator (e.g., Barklem et al. 
2002). Adopting the theory of Ali & Griem 
(1965, 1966) for resonance broadening and 
Griem (1960) for Stark broadening, we 
computed Ha profiles for several effective 
temperatures, using the code SYNTHE 
(Kurucz et al. 2005; Sbordone et al. 2005). 
For further details on the computations of 
hydrogen lines in SYNTHE, see Castelli 
& Kurucz (2001) and Cowley & Castelli 
(2002). 

Figure [2] compares these synthetic Ha 
profiles with the observed profile for sev- 
eral temperatures. To evaluate the good- 
ness of fit, we employ the same method 
as in the previous section. However, in 
this case M = 1 (the effective tempera- 
ture T e g), fi is the synthetic normalized 
flux, Fi is the observed normalized flux, 
and (jj = 1/SNR. The SNR was estimated 
constant average value in continuum 
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Fig. 2. — Upper panel: Synthetic Ha profiles for three effective temperatures compared to 
the observed spectrum of Tycho G (SNR « 30) and normalized to the level of the observed 
spectrum of Tycho G at 6540 A. We also display the fitting regions at the top. Lower panel: 
Results of x%u analysis for three continuum positions separated by ±1/SNR = 0.03: 1.03 
(left), 1.00 (middle), and 0.97 (right). 



regions close to the observed Ha profile. 
The fitting regions are indicated in Fig- 
ure [21 which contains all of the spectral 
regions close the center of the Ha pro- 
file where there are no stellar lines and 
where the normalized flux is greater than 
~ 0.7. The best fit provided a temperature 
of T cfr = 6050 ± 250 K, with the error bars 
estimated from the best fits after displac- 
ing the observed continuum up and down 
by 1/SNR = 0.03. 

In the following section, we adopt the 
currently most reliable method for the de- 
termination of the stellar parameters. 



3.3. Ionization Equilibrium of Iron 

The best determination of the atmo- 
spheric parameters of Tycho G can be ob- 
tained from the ionization equilibrium of 
iron. We measured the equivalent widths 
of 32 Fe I and 10 Fe II isolated lines using 
an automatic line-fitting procedure which 
performs both line detection and Gaus- 
sian fits to unblended lines (Frangois et 
al. 2003). The analysis was done us- 
ing the 2002 version of the code MOOG 
(Sneden 1973) and a grid of local ther- 
modynamic equilibrium (LTE) model at- 
mospheres (Kurucz 1993). We adopted 
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Fig. 3. — Fe I and Fe II abundances computed using the best-fit atmospheric parameters of 
Tycho G (left panels) and two sets of stellar parameters: T c g = 5900 K and log(g/cm s2) 
= 2 dex (middle panels), T e g = 6200 K and log((?/cm s2) = 4.35 dex (right panels). Upper 
panels: Fe I abundance vs. excitation potential; middle panels: Fe I abundance vs. reduced 
equivalent widths; bottom panels: Fe II abundance vs. excitation potential. The dashed 
lines in the upper and middle panels represent the fits to the points, and in the bottom 
panels the mean abundance. 



the atomic data from Santos, Israelian, 
& Mayor (2004), where the log gf values 
are adjusted until the solar atlas (Kurucz 
et al. 1984) is reproduced with a Kurucz 
model for the Sun having T e g- = 5777 K, 
log(#/cm s2) = 4.44 dex, & = 1.00 km s _1 , 
and log e(Fe) = 7A7 dex. This line list 
was designed to determine stellar parame- 
ters of planet-host stars with roughly the 
same atmospheric parameters and metal- 
licity as Tycho G; it provides a compila- 



tion of 42 completely unblended Fe lines 
ideal for the determination of stellar atmo- 
spheric parameters. Note that none of the 
stellar lines used in this work for chemical 
analysis is affected by telluric lines. 

The atmospheric parameters were de- 
termined by iterating until correlation co- 
efficients between loge(Fe I) and xu as 
well as between loge(Fe I) and \og{W\/\), 
were zero, and the mean abundances from 
Fe I and Fe II lines were the same (see left 
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panels of Fig. [3]) . The derived parameters 
were T eff = 5900 ± 100 K, log(#/cm s2) = 
3.85±0.30 dex, & = 1.23±0.23 km s _1 , and 
[Fe/H] = -0.05 ± 0.09. The uncertainties 
in the stellar parameters were estimated as 
described by Gonzalez & Vanture (1998). 
Thus, the uncertainties in £ t and T eS take 
into account the standard deviation of the 
slope of the least-squares fits, loge(Fe I) 
versus \og(W\/\) and loge(Fe I) versus xi- 
The uncertainty in log g considers the 
uncertainty on T e g in addition to the scat- 
ter of the Fe II abundances. Finally, the 
uncertainty in Fe abundance is estimated 
from a combination of the uncertainties in 
£t and T eff , in addition to the scatter of the 
Fe I abundances, all added in quadrature. 

In Figure [31 we also show two set of 
stellar parameters, T eS = 5900 K and 
log (g /cm s2) = 2 dex (middle panels), and 
T eff = 6200 K and log(#/cm s2) = 4.35 dex 
(right panels), corresponding to the spec- 
tral types F8 I and F6-7 V (Gray 1992). 
The lower-middle and upper-middle panels 
show inconsistent abundances between Fe I 
and Fe II, whereas the upper-right panel 
shows a negative slope for iron abundance 
versus excitation potential xi °f the spec- 
tral lines Fe I. Therefore, these two sets of 
parameters cannot be the stellar parame- 
ters of Tycho G. According to the spectral 
classification provided by Gray (1992), the 
stellar parameters of Tycho G correspond 
to a spectral type of G0-1 IV. 

Non-LTE (NLTE) corrections for stars 
with similar stellar parameters and metal- 
licity have been estimated at < 0.07 dex 
for Fe I lines, while Fe II seems to be insen- 
sitive to NLTE effects (Thevenin & Idiart 
1999). In addition, Santos et al. (2004) 
compared spectroscopic surface gravities 



with log g using Hipparcos parallaxes and 
found an average difference of ~ +0.03 
dex, which might in fact reflect NLTE ef- 
fects on Fe I lines. Our assumption of LTE 
seems to be accurate enough to determine 
the surface gravity of the star within the 
error bars. Note that ionization equilib- 
rium also holds for Si and Cr (see § HJ). 

This determination of effective tempera- 
ture is consistent within the error bars with 
that derived from the Ha profile (T c g = 
6050 ± 250 K). However, echelle spectro- 
graphs are not ideally suited for measur- 
ing the shape of broad features such as the 
wings of Balmer lines (Allende Prieto et 
al. 2004), and the SNR of ~30 is not high 
enough for the determination of effective 
temperature from the Ha profile with great 
accuracy. For this reason, we will adopt 
the T cfT = 5900 ± 100 K value for the de- 
termination of the chemical abundances. 

4. Chemical Abundances 

Using the derived stellar parameters, we 
determine the element abundances from 
the absorption-line equivalent width (EW) 
measurements of the elements given in Ta- 
ble HI except for lithium and sulfur for 
which we use a %2 procedure to find the 
best fit to the observed features (e.g., 
Gonzalez Hernandez et al. 2004, 2005, 
2006). The atomic data for O, C, S, and 
Zn are taken from Ecuvillon, Israelian, & 
Santos (2004, 2006), those of Na, Mg, Al, 
Si, Ca, Sc, Ti, V, Cr, Mn, Co, and Ni are 
from Gilli, Israelian, & Ecuvillon (2006), 
and those for Sr, Y, and Ba are from Reddy 
et al. (2003). 

In Figures HI [5], and [6] we display sev- 
eral spectral regions of the observed spec- 
trum of Tycho G in comparison with syn- 
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Fig. 4. — Best synthetic spectral fits (solid lines) to the Keck HIRES spectrum of Tycho G 
for two spectral regions. The core of the Ha line suffers from NLTE effects that are not 
taken into account in our LTE synthetic spectrum. 



thetic spectra computed using the derived 
element abundances. Uncertainties in the 
abundances of all elements were then de- 
termined, adding in quadrature the er- 
rors due to the sensitivities of the result- 
ing abundances to changes in assumed at- 
mospheric parameters and the dispersion 
of the abundances from individual lines of 
each element. 

In Table [T] we provide the average abun- 
dance of each element together with the er- 
rors. The errors in the element abundances 
show their sensitivity to the uncertainties 
in the effective temperature (A Teff ), surface 



gravity (Ai ogs ), microturbulence (A^), and 
the dispersion of the measurements from 
different spectral features (A CT ). The er- 
rors A CT were estimated as A CT = cr/v^/V, 
where a is the standard deviation of the N 
measurements. 

The abundances of several elements 
listed in Table [I] relative to iron are com- 
pared in Figures ElIEJ and El with the Galac- 
tic trends of these elements in the relevant 
range of metallicities. The abundances of 
some heavy elements relative to iron are 
above the solar values. The value found 
for [Ni/Fe], 0.16 ± 0.04, is well above the 
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Table 1 

Chemical Abundances of Tycho G 



Species loge(X) a [X/H] [X/Fc] a A CT A Trff A log9 A e A[X/H] A[X/Fc] n b 



C I 


8.56 


-0.15 


-0.07 


0.15 


0.07 


-0.06 


0.09 





0.13 


0.19 


5 


O I 


8.74 


-0.02 


0.03 


0.02 


0.01 


-0.08 


0.07 


-0.02 


0.11 


0.18 


3 


Na I 


6.33 


0.04 


0.09 


0.18 


0.10 


0.05 


-0.05 


-0.02 


0.12 


0.12 


3 


Mgl 


7.58 


0.08 


0.13 


0.06 


0.03 


0.05 


-0.04 


-0.03 


0.08 


0.06 


4 


Al I 


6.47 


0.16 


0.21 


0.05 


0.03 


0.06 


-0.01 


-0.01 


0.07 


0.05 


2 


Si I 


7.55 


-0.03 


0.02 


0.09 


0.02 


0.04 


-0.01 


-0.02 


0.05 


0.05 


12 


Si II 


7.55 


-0.03 


0.02 






-0.08 


0.11 


-0.03 


0.14 


0.20 


6371 


S I c 


7.21 


-0.02 


0.03 


0.07 


0.05 


-0.03 


0.12 





0.13 


0.18 


2 


Ca I 


6.36 


-0.04 


0.01 


0.11 


0.03 


0.08 


-0.04 


-0.05 


0.11 


0.05 


10 


Sc II 


3.10 


-0.06 


-0.01 


0.09 


0.03 


0.01 


0.12 


-0.04 


0.13 


0.15 


7 


Ti I 


4.99 


-0.03 


0.02 


0.14 


0.04 


0.11 





-0.02 


0.12 


0.06 


10 


V I 


4.00 


0.07 


0.12 


0.21 


0.07 


0.10 


-0.01 


-0.01 


0.12 


0.08 


8 


Cr I 


5.67 


-0.04 


0.01 


0.16 


0.05 


0.07 





-0.01 


0.09 


0.06 


9 


Cr II 


5.67 


-0.03 


0.02 






-0.02 


0.11 


-0.04 


0.12 


0.16 


5305 


Mn I 


5.39 


0.02 


0.07 


0.15 


0.05 


0.09 


-0.02 


-0.05 


0.12 


0.06 


8 


Fe I 


7.47 


-0.05 





0.12 


0.02 


0.08 


-0.01 


-0.04 


0.09 




32 


Fe II 


7.47 


-0.05 





0.10 


0.03 





0.13 


-0.05 


0.14 


0.17 


10 


Co I 


4.92 


0.13 


0.18 


0.14 


0.07 


0.09 





-0.02 


0.12 


0.08 


4 


Ni I 


6.25 


0.11 


0.16 


0.19 


0.04 


0.08 


-0.01 


-0.04 


0.10 


0.04 


22 


Zn I 


4.60 


-0.06 


-0.01 


0.09 


0.05 


0.03 


0.02 


-0.08 


0.10 


0.09 


3 


Sr I 


2.64 


0.07 


0.12 






0.09 


-0.02 


-0.06 


0.11 


0.03 


4607 


Y II 


2.12 


-0.14 


-0.09 


0.16 


0.08 


0.02 


0.12 


-0.06 


0.16 


0.17 


4 


Ba II 


2.20 


0.19 


0.24 


0.12 


0.07 


0.03 


0.05 


-0.17 


0.19 


0.17 


3 



Note. — Chemical abundances of Tycho G and uncertainties produced by A(T e g) = +100 K, A(logg) = +0.3 
dcx, and A(£) = +0.23 kms -1 . 

a The solar element abundances were adopted from Santos ct al. (2004), Ecuvillon ct al. (2004, 2006), Gilli et 
al. (2006), and Reddy et al. (2003). 

b Number of spectral lines of this element analyzed in the star, or if there is only one, its wavelength. 

c These abundances were determined by fitting the observed spectra with synthetic spectra computed with the 
LTE code MOOG. 
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Fig. 5. — Same as Fig. HI but for other spectral regions. 



average value if we consider the intrinsic 
dispersion ([Ni/Fe] = -0.07 ± 0.04, see 
Fig. [9]) for the stars with similar metallic- 
ity of Tycho G ([Fe/H] = -0.05 ± 0.09). 
In Figure [5] we display several spectral 
ranges where some Ni I lines are well repro- 
duced. Generally, Ni tracks Fe throughout 
the [Fe/H] range down to [Fe/H] = —1. 
The average value and scatter found by 
Reddy, Lambert, & Allende Prieto (2006) 
is [Ni/Fe] = -0.05 ± 0.02 for thin-disk 
stars and [Ni/Fe] = -0.01 ±0.04 for thick- 
disk stars. Bensby et al. (2005) found 
[Ni/Fe] = -0.06 ± 0.04 for thin-disk stars 
and [Ni/Fe] = -0.02 ± 0.02 for thick-disk 
stars. This suggests pollution in Ni by Ty- 



cho SN 1572 ejecta. We find in general 
solar or slightly above solar abundances of 
heavy elements, while the a elements are 
consistent with the solar value or below. 

The binary companions of black holes 
or neutron stars show a different enhance- 
ment of metals with respect to the so- 
lar values. These stars exhibit signifi- 
cant enrichment in a elements. For in- 
stance, the black hole binary Nova Scor- 
pii 1994 ([Fe/H] = -0.1) shows enhance- 
ments of [a/Fe] = 4-8 in Mg, S, Si, and O 
(Gonzalez Hernandez et al. 2008a). The 
pollution seems to be related to the fall- 
back of the ejected material of the super- 
nova. In SNe la, one does not expect fall- 
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Fig. 6. — Same as Fig. HI but for other spectral regions. 



back of material onto the companion star 
since the compact object (the white dwarf 
in this case) is destroyed and the gravita- 
tional potential well is not sufficiently deep 
to retain the ejected material. Therefore, 
one expects to have a low contamination in 
intermediate-mass elements. The material 
captured by the companion should consist 
of heavy elements; they are less likely to es- 
cape the companion star since they move 
at lower velocities in the supernova ejecta. 
A discussion of this is given in § 18. 1[ 

Another feature of interest in Tycho G 
is its high lithium abundance; the Li line 
at 6708 A is pronounced. Modeling the 
Li abundance gives A(Li) = 2.50 ± 0.09. 



The lithium abundance is provided here 
as A (Li) = log[iV(Li)/jV(H)] + 12 and 
takes into account NLTE effects. The 
NLTE abundance correction, Anlte = 
loge(X) NLTE - loge(X) LTE , for Li was de- 
rived from the theoretical LTE and NLTE 
curves of growth in Pavlenko & Magazzu 
(1996). 

The high Li abundance in Tycho G 
is intriguing since Li is easily destroyed 
in the convective envelopes of stars that 
have evolved away from the main se- 
quence. Subgiants with A(Li) > 2.2 are 
rare. Thoren, Edvardsson, & Gustafsson 
(2004) found that subgiants with detected 
Li divide into two groups: the stars with 
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Fig. 7. — Abundance ratios of Tycho G (wide cross) in comparison with the abundances of 
G and K metal-rich dwarf stars. Galactic trends were taken from Ecuvillon et al. (2004), 
Ecuvillon et al. (2006), and Gilli et al. (2006). The size of the cross indicates the la 
uncertainty. For the abundance of oxygen in metal-rich dwarfs, we have only considered 
abundance measurements in NLTE for the triplet O I 7771-5 A. The dashed-dotted lines 
indicate solar abundance values. 



T eff < 5500 K have A(Li) < 1.5, and a 
warmer group of stars have higher abun- 
dances, although still much lower than 
their primordial abundances, with A (Li) 
< 2.2. This seems to go in the direction of 
confirming post-main-sequence evolution 
toward lower amounts of Li as expected 
from stellar models. 

Tycho G shares a high Li content with 
a known sample of binary companions of 
neutron stars and black holes (Martin et 
al. 1994a; Gonzalez Hernandez et al. 2004, 



2005). However, we know that cataclysmic 
variables, whose compact objects are white 
dwarfs, do not show such high Li abun- 
dances (Martin et al. 1995). We discuss 
this further in § 18.21 

5. Distance 

One can estimate the distance of Ty- 
cho G from different photometric colors 
and stellar parameters. We use the photo- 
metric magnitudes in five different filters: 
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Fig. 8. — Same as Fig. [3, but for different elements. 



m v = 18.71 ± 0.04 and m R = 17.83 ± 0.03 
mag from Ruiz-Lapuente et al. (2004), and 
mj = 15.84±0.07, m H = 15.16±0.11, and 
m K = 15.03 ± 0.12 mag from the 2MASS 
catalog^]. We derived the radius of the 
star from the surface gravity, log((7/cm s2) 
= 3.85 ± 0.3 dex, and assuming a mass 
of 1 M . This radius, together with the 
spectroscopic estimate of the effective tem- 
perature, T eff = 5900 ± 100 K, provides 
an intrinsic bolometric luminosity of 1.9 < 



lr The Two Micron All Sky Survey is a joint project 
of the University of Massachusetts and the In- 
frared Processing and Analysis Center/California 
Institute of Technology, funded by the National 
Aeronautics and Space Administration (NASA) 
and the National Science Foundation (NSF). 



LJL® < 7.6. 

In Table El we show distance deter- 
minations for the different filters, using 
the bolometric corrections for models with- 
out overshooting (Bessell, Castelli, & Plez 
1998), for our preferred value of the color 
excess, E(B — V) = 0.76 mag, and other 
different sets of the relevant parameters. 
We compute the magnitude corrected for 
extinction in each filter as m^o = my— Ay, 
where Ay = 3.12 E(B — V) is the extinc- 
tion in the Johnson V filter. We adopt the 
following values for other filters (Schaifers 
et al. 1982): A R /A V = 0.84, Aj/A v = 
0.32, A H /A V = 0.27, and A K /A V = 0.21. 

Table [2] clearly shows the uncertainty 
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Fig. 9. — Same as Fig. [3, but for different elements. 



in the distance due to the uncertainty in 
the stellar parameters of Tycho G. Note 
that the relatively small formal error in 
the distance estimate for each filter is cal- 
culated by assuming that the magnitudes 
equal rrii + Am;. 

In Table [2j we also show for each set 
of parameters the average distance de- 
terminations weighted by the uncertainty 
of each individual distance determination. 
These values cover a distance range 2.48- 
4.95 kpc with errors of 0.21-0.52 kpc. 
We have also considered the possibility 
that Tycho G lost part of its envelope, 
AM = -0.2 M , due to the impact of 
the supernova shock wave (Marietta, Bur- 
rows, & Fryxell 2000; Pakmor et al. 2008). 



However, a lower mass does not produce 
a significant change in the derived aver- 
age distance. In all cases, we have as- 
sumed that the companion star is able to 
almost completely recover thermal equi- 
librium between ~ 102 and ~ 103 yr 
after the white dwarf explosion (Podsi- 
adlowski 2003). The distance to the SN 
1572 remnant inferred from the expansion 
of the radio shell and by other methods 
is 2.83 ± 0.79 kpc (Ruiz-Lapuente 2004). 
Such a distance is in good agreement with 
the derived distances of Tycho G. 
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Table 2 
Distance of Tycho G 



Parameter 


E(B - V) a = 0.76 E(S - 


V) & = 0.66 T cS b 


= 5800 K 


log g c = 4.15 


log g d = 3.55 


M* e = 0.8 M 


il*(R ) 


1.87 


1.87 


1.87 


1.32 


2.63 


1.67 


L* (L Q ) 


3.80 


3.80 


3.54f 


1.90 


7.57 


3.04 


d mv , f (kpc) 


3.80 ±0.07 


4.39 ± 0.08 


3.65 ± 0.07 


2.69 ±0.05 


5.38 ± 0.10 


3.40 ±0.06 


d m R (kpc) 


3.52 ±0.05 


3.98 ± 0.06 


3.41 ± 0.05 


2.49 ± 0.03 


4.99 ± 0.07 


3.15 ± 0.04 


dm j (kpc) 


3.49 ±0.11 


3.65 ± 0.12 


3.43 ± 0.11 


2.47 ±0.08 


4.94 ±0.16 


3.12 ± 0.10 


dm H (kpc) 


3.09 ±0.15 


3.22 ± 0.16 


3.06 ± 0.15 


2.20 ±0.11 


4.37 ± 0.22 


2.77 ±0.14 


(kpc) 


3.16 ±0.18 


3.25 ± 0.18 


3.13 ± 0.18 


2.24 ±0.13 


4.46 ±0.25 


2.82 ± 0.16 


dl v (kpc) 


3.50 ±0.31 


3.85 ± 0.52 


3.40 ±0.25 


2.48 ±0.21 


4.95 ± 0.44 


3.13 ± 0.27 


a The distance 


is estimated by adopting T c g 


= 5900 K, log(g/cm 


s2) = 3.85 dcx, M* = 1 Mq, and the color excess 


indicated in each 



column by E(B ~ V). 

b The distance is estimated by adopting T eS = 5800 K, log(#/cm s2) = 3.85 dcx, M* = 1 Mq, and E(B - V) = 0.76 mag. 
c The distance is estimated by adopting T eff = 5900 K, log(#/cm s2) = 4.15 dex, M* = 1 M , and E(B - V) = 0.76 mag. 
d The distance is estimated by adopting T eff = 5900 K, log(#/cm s2) = 3.55 dex, M* = 1 M , and E(B - V) = 0.76 mag. 
e Thc distance is estimated by adopting T cff = 5900 K, log(#/cm s2) = 3.85 dcx, M* = 0.8 M , and E(B - V) = 0.76 mag. 
'The error bar in the distance is estimated from the uncertainty in the magnitude. 

g Average distance weighted by the uncertainty of each individual distance determination. The error bar shows the large dispersion among 
the individual determinations. 



6. Projected Rotational Velocity of 
Tycho G 

Tycho G is not a fast rotator, unlike the 
companion stars of Type II SNe. The pro- 
jected rotational velocity is constrained by 
the instrumental resolution to v ro t sin i < 
6.6 km s -1 . 

However, the companion stars of SNe II 
and other core-collapse SNe have different 
rotational histories from those proposed for 
SN la progenitors; they would be massive 
stars, and maybe fast rotators (Meynet 
& Maeder 2005, and references therein). 
Moreover, the orbital elements just before 
the explosion are also different; in par- 
ticular, the separations are larger. Given 
the orbital history up to the explosion, it 
is reasonable to find fast rotators around 
black holes and neutron stars (Gonzalez 
Hernandez et al. 2004, 2005, 2008b). The 



companion stars of these compact objects 
do not transfer mass to the star that be- 
comes a core-collapse supernova. If the ex- 
plosion does not unbind the system, the 
companion can be seen as a star orbiting 
a black hole or neutron star and show- 
ing significant rotational velocity. But, for 
SNe la, we should not expect such high 
velocities since the companion of the white 
dwarf would have transferred mass and an- 
gular momentum to the white dwarf. And, 
after the impact, the loss of the envelope 
would have further helped decrease the an- 
gular momentum and to brake rotation. 

Let us consider, for instance, a WD 
mass Mwd = 1-4 M and a companion 
of mass M = 1 M Q (a system similar to 
U Scorpii; see, e.g., Thoroughgood et al. 
2001). In such a system, an orbital veloc- 
ity of 90 km s -1 (this value could roughly 
correspond to the peculiar velocity of Ty- 
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cho G estimated by Ruiz-Lapuente et al. 
2004, still consistent with the values from 
the present work) corresponds to an orbital 
separation a = 19.28 R and an orbital 
period P orb = 6.32 d. The effective Roche 
lobe radius of the companion can be cal- 
culated from the Eggleton (1983) approx- 
imation (accurate up to a few percent for 
all mass ratios): 



L [0.6 + g- 2 / 3 ln(l + gV3)J ' v ; 

It would then be R L = 6.75 R . 

Assuming that tidal interaction had 
been effective enough to keep the rota- 
tional period of the companion locked with 
the orbital period in spite of the angu- 
lar momentum loss due to mass transfer, 
and that no further changes have inter- 
vened, the rotation period of the compan- 
ion should now be P TOt = P OI h = 6.32 d. 
Therefore, the rotational velocity is v TOt = 
2TcR/P rot . However, the present radius of 
the Tycho G is within the range R ~ 1-2 
R , as can be seen from its position in the 
Hertzsprung-Russell diagram. This then 
gives f ro t ~ 8-16 kms -1 . Here we have 
assumed that the impact of the SN ejecta 
has removed the outer, more distended and 
less gravitationally bound layers (Marietta 
et al. 2000; Pakmor et al. 2008), and that 
the star is able to completely or partially 
recover thermal equilibrium ~ 500 yr af- 
ter the explosion (Podsiadlowski 2003). In 
this case, when the energy deposited on 
the star due to the impact of the shock 
wave of the supernova is < 5 x 10 46 erg, 
the luminosity of the star in the time inter- 
val ~ 102-103 yr after the explosion could 
be roughly ~ 2-5 L . Hence, the radius 



of the star would be in the range ~ 1.5-2 
R . 

Thus, an observed value of v rot sin i < 
6.6 km s^ 1 is not an unexpected result in 
this case (for a typical angle % « 45° or 
lower). In fact, the preceding analysis does 
not even take into account angular momen- 
tum losses due to mass stripping, nor the 
slowing down of rotation (with conserva- 
tion of angular momentum) that would re- 
sult from puffing up of the fraction of the 
envelope remaining bound. 

Any suggestion that the rotational ve- 
locity of the former companion of a SN la 
should now be very high is based on the 
possibly incorrect assumption that its cur- 
rent radius should remain, after the SN la 
explosion, at the large value it had at the 
time of Roche-lobe filling (equal to 6 R in 
our example). Thus, expected rotational 
velocities of companions of supernovae, as 
well as potential pollution by the super- 
nova ejecta, are different in SNe la and 
SNe II/Ib/Ic. 

7. Space Velocity and Position of 
Tycho G 

In our previous study (Ruiz-Lapuente et 
al. 2004), we measured the proper motion 
of Tycho G to be fib = —6.11 ± 1.34 mas 
yr -1 and fii = —2.6 ± 1.34 mas yr -1 . The 
proper motion in Galactic latitude implies 
that the star is an outlier in proper motion, 
with a derived tangential velocity of 94±27 
kms" 1 . A new, more accurate and precise 
proper motion is currently being derived 
from HST images obtained in Cycle 16. 

The spectral data presented by Ruiz- 
Lapuente et al. (2004) allowed an es- 
timate of the heliocentric radial velocity 
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of Tycho G with an uncertainty of 10-18 
kms" 1 . The new Keck HIRES data pro- 
vide a much smaller uncertainty. For this 
new measurement we selected six strong, 
unblended lines in the spectral range 6000- 
8000 A and measured the radial velocity 
of each individual line. The average radial 
velocities from those lines are found to be 
-87.36 ± 0.60 kms" 1 and -87.53 ± 0.48 
km s _1 during the first and second nights, 
respectively, where the errors indicate the 
dispersion of the measurements. The aver- 
age of the two epochs is v r = —87.4 ± 0.5 
kms' 1 . This value implies that v r ~ 
—80 km s _1 in the local standard of rest 
(LSR). In our previous measurement, we 
had quoted v r = —99 ± 6 km s" 1 (statis- 
tical error). The measurement was made 
using Ha and H/5 lines, which are worse 
indicators of the velocity of the star. 

Stars in the direction of the Tycho SN 
remnant at distances of 2-4 kpc move at 
an average radial velocity v r ~ —20 to 
—40 kms" 1 (in the LSR), in agreement 
with Galactic rotation models (Brand & 
Blitz 1993; Dehnen & Binney 1998). Ty- 
cho G is well above the average of stars 
in its vicinity. Both the radial velocity and 
the proper motion show a comparably high 
velocity. One could argue that Tycho G 
might be at a larger distance (although this 
seems unlikely from its photometric mag- 
nitudes; see §[5]), but the observed proper 
motion would make the tangential velocity 
unreasonably high if the star were > 6 kpc 
away (vt ~ 187 km s" 1 at d — 6 kpc). On 
the other hand, a distance of 4-5 kpc (al- 
though v t ~ 140 km s" 1 at d — 4.5 kpc) 
might not be ruled out by these arguments, 
albeit somewhat inconsistent with the es- 
timated distance of the Tycho SN 1572 
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Fig. 10. — Upper panel: Toomre diagram 
of the thick-disk stars (black filled circles), 
thin-disk stars (open circles), and transi- 
tion stars (grey filled circles) from Bensby 
et al. (2003, 2005). The star Tycho G is 
displayed as a rhombus. Lower panel: the 
same as upper panel, but for stars with 
metallicity greater than or equal to the 
metallicity of Tycho G. 

remnant (2.83 ± 0.79 kpc; Ruiz-Lapuente 
2004). 

However, Fuhrmann (2005) has ar- 
gued, based on the Galactic space ve- 
locity components of Tycho G relative 
to the LSR ({7 L sr, Vlsr, Wlsr), that this 
star could be just a thick-disk star pass- 
ing close to the SN remnant. Assum- 
ing that the Sun moves with components 
(^0,lsr,K 3 , LS r,W ,lsr) = (10., 5.2, 7.2) 
km s _1 (Dehnen & Binney 1998) relative to 
the LSR, for distance of 3 kpc, we estimate 
the space Galactic velocity components 
of Tycho G to be (U LSK ,V LSR ,W LSK ) = 
(85 ± 17, -50 ± 10, -82 ± 24) km s _1 . In 
Figure [10] we display Galactic velocity of 
Tycho G in comparison with the Galactic 
velocities of thick-disk and thin-disk stars. 
Although the Galactic velocity of Tycho G 
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is comparable to that of thick-disk stars, 
its metallicity is more common in thin-disk 
stars. In the lower panel of Figure [TOj we 
display only those stars with metallicity 
equal to or greater than the metallicity of 
Tycho G. The number of thick-disk stars is 
reduced considerably, but there are still a 
few thick-disk stars with the same metal- 
licity and Galactic velocity as Tycho G. On 
the other hand, the high Ni abundance of 
Tycho G would be very unusual for a thick- 
disk star. We therefore conclude that while 
a thick-disk origin may be consistent with 
the data, this possibility is not favored be- 
cause it provides no explanation for the 
anomalous Ni abundance. 

Although the position of Tycho G in the 
sky does not coincide with the proposed 
geometrical center of the X-ray and radio 
remnants of SN 1572, the offset of the po- 
sition is within the uncertainties due to 
the asymmetry resulting from interaction 
of the expanding remnant with the circum- 
stellar medium. Indeed, there is 0.56' dis- 
placement along the E-W axis between the 
radio emission and the high-energy contin- 
uum observed by XMM- Newton, in the po- 
sition of the western rim (Decourchelle et 
al. 2001). Such asymmetry amounts to a 
14% offset along the E-W axis. It appears 
to be due to the ejecta having encountered 
a dense H cloud at the eastern edge, giv- 
ing rise to brighter emission and lower ex- 
pansion velocity of the ejecta there, while 
interacting with a lower-density medium 
in the western rim. It was precisely ow- 
ing to those considerations that we stud- 
ied all of the 23 stars with my < 22 
mag within a circle of radius 0.65' around 
the geometrical center of the X-ray emis- 
sion as given by the Chandra observations 



(Hughes 2000). All of the stars within that 
circle were considered on an equal footing 
as possible companions until all were dis- 
carded with the exception of Tycho G, due 
to the unique characteristics of this star. In 
the next two sections, we discuss its chem- 
ical peculiarities. 

8. Discussion 

8.1. Has Tycho G been Polluted by 
the SN la? 

The measured ratio [Ni/Fe] = 0.16 ± 
0.04 in Tycho G is 4-5cr above the average 
value of this ratio in the Galaxy ([Ni/Fe] 
~ —0.05). It indicates an overabundance 
of Ni with respect to the Galactic trend. 
Ni may have been trapped by the star as 
slowly moving material from deep layers of 
the SN ejecta orbited around it. While the 
rapidly moving ejecta consisting mainly of 
intermediate-mass elements would have es- 
caped away from the site of the explosion, 
the very slowly moving tail made of heavy 
elements like Ni could have been captured 
by the companion. 

To be consistent with the evolutionary 
path suggested above (and also by Ruiz- 
Lapuente et al. 2004), we use a = 19.28 R 
for the separation between the white dwarf 
and its companion, and R* = Rl = 6.75 
R for the radius of the companion (just 
before the SN explosion). 

The mass that could be trapped by the 
companion is given by the subtended an- 
gle, 

m t = AM (ttRI/Att a 2 ) /, (3) 

where AM is the ejected mass and / is the 
fraction of retained material. 



19 



Table 3 

Supernova Ia models in Tycho G 



[X/H] EXPECTED 



ELEMENT 


[X/H] OBSERVED a 


[X/H] » 


W7 


W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


O 


-0.02 ±0.11 


-0.15 


-0.14 


-0.14 


-0.14 


-0.15 


-0.15 


-0.14 


-0.15 


Na 


0.04 ±0.12 


-0.26 


-0.26 


-0.26 


-0.26 


-0.26 


-0.26 


-0.26 


-0.26 


Mg 


0.08 ±0.08 


-0.10 


0.04 


0.13 


0.03 


-0.02 


-0.05 


0.03 


-0.02 


Al 


0.16 ±0.07 


-0.06 


-0.05 


-0.06 


-0.06 


-0.06 


-0.06 


-0.06 


-0.06 


Si 


-0.03 ±0.05 


-0.22 


-0.07 


-0.08 


0.02 


-0.03 


-0.06 


0.02 


-0.03 


s 


-0.02 ±0.13 


-0.26 


-0.10 


-0.10 


-0.00 


-0.05 


-0.09 


0.00 


-0.06 


Ca 


-0.04 ±0.11 


-0.31 


-0.15 


-0.08 


0.03 


-0.03 


-0.08 


0.03 


-0.03 


Ti 


-0.03 ±0.12 


-0.06 


0.00 


0.02 


0.12 


0.10 


0.09 


0.08 


0.07 


Cr 


-0.04 ±0.09 


-0.27 


0.03 


0.05 


0.27 


0.21 


0.16 


0.24 


0.20 


Mn 


0.02 ±0.12 


-0.27 


0.12 


0.05 


0.11 


0.06 


0.03 


0.09 


0.05 


Fe 


-0.05 ±0.09 


-0.21 


0.11 


0.12 


0.09 


0.13 


0.15 


0.08 


0.14 


Co 


0.13 ±0.12 


-0.14 


0.03 


0.02 


-0.07 


-0.03 


-0.01 


-0.09 


-0.03 


Ni 


0.11 ±0.10 


-0.27 


0.40 


0.32 


0.06 


0.16 


0.23 


0.04 


0.18 



a Observed abundances of the secondary star in Tycho G. 

b Initial abundances assumed for the secondary star in Tycho G; see text. 

c Expccted abundances in Tycho G after being contaminated with nucleosynthetic products of SN Ia explosion models from 
Iwamoto et al. (1999). Some of them are presented in Fig. 1111 



We try to model the contamination of 
the companion star by the nucleosynthetic 
products of different SN Ia models from 
Iwamoto et al. (1999). In Table [3] we show 
the expected abundances of Tycho G af- 
ter having captured a significant amount 
of the SN ejecta. In these model compu- 
tations, we adopt the yields from SN Ia 
models with different kinetic explosion en- 
ergies in the range (1.30-1.44) x 10 51 ergs, 
different central densities (in units of 109 
g cm" 3 ) of p 9 = 1.37 (C) and 2.12 (W), and 
different deflagration velocities with fast 
deflagration in the models W7 and W70 
(Nomoto, Thielemann, & Yokoi 1984) and 
slow deflagration in the others. We as- 
sumed as initial abundances the average 
abundances of disk stars with a metallic- 
ity of [Fe/H] = -0.20 ± 0.09, and a cap- 
ture efficiency factor / = 0.03. The mat- 
ter that is captured by the companion has 



a much larger mean molecular weight than 
the composition of its atmosphere and then 
is completely mixed with the whole star 
due to thermohaline mixing, in a short 
timescale (Podsiadlowski et al. 2002). 

Badenes et al. (2006) found, by com- 
paring a grid of X-ray synthetic spectra 
based on hydrodynamical models, that 
the fundamental properties of X-ray emis- 
sion in Tycho SN 1572 are well repro- 
duced with one-dimensional delayed det- 
onation models having a kinetic energy of 
~ 1.2 x 10 51 ergs. In Table El the mod- 
els called WDD1,2,3 and CDD1,2 include 
transitions from slow deflagration to det- 
onation, which are equivalent to delayed 
detonation models. 

In Figure [TU we display the observed 
abundances of Tycho G in comparison 
with the expected abundances from sev- 
eral SN Ia models. Except for Na and 
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Al, we find reasonable agreement between 
the observed abundances of Tycho G and 
the expected abundances after contami- 
nation from the nucleosynthetic products 
of the SN la. Surprinsingly, Gonzalez 
Hernandez et al. (2008a) measured the 
chemical abundances of the secondary star 
of the black hole binary Nova Scorpii 1994 
and found several a elements significantly 
enhanced in the companion star while Al 
and Na were not enhanced. The enhanced 
a elements in that star were explained 
by comparing the observed abundances 
with supernova/hypernova explosion mod- 
els. However, these models also predict 
that Al and Na should be enhanced at a 
comparable level as oxygen, and they did 
not find an explanation for this discrep- 
ancy. Therefore, the model of contamina- 
tion of Tycho G by using SN la models 
could provide an explanation for the rela- 
tively high observed [Ni/Fe] with respect 
to the Galactic trend of this element. 



0.4 



0.2 



0.0 



-0.2 



-0.4 




O NaMeAl Si S 



10 



Ca Ti CrMnFeCoNi 
20 25 



Fig. 11. — Expected abundances in Tycho 
G after contamination from the nucleosyn- 
thetic products in SN la models (Iwamoto 
et al. 1999), in comparison with the ob- 
served abundances. 



One could find stable isotopes of Co as 
well among the slowly moving material in 
a SN la. The ratio [Co/Fe] = 0.18 ± 0.08 
in Tycho G is consistent with contamina- 
tion from the SN ejecta. The study by 
Reddy et al. (2006) finds that [Co/Fe] 
shows flat behavior with [Fe/H] for thin- 
disk stars ([Co/Fe] = -0.05 ± 0.02). For 
thick-disk stars, a weak trend of decreas- 
ing [Co/Fe] with increasing [Fe/H] is ap- 
parent, falling to [Co/Fe] = 0.00 ± 0.04 at 
the metallicity of Tycho G. The effect of 
a larger [Co/Fe] is > la, a lower signifi- 
cance level compared to the ratio [Ni/Fe] 
(see also Fig. 6). The ratios of the a el- 
ements (O, Si, Ca, and Ti) to iron do not 
show enhancements over Galactic values. 

A study of the SN 1572 remnant in X- 
rays (Hamilton et al. 1985) suggests that it 
contains ~0.66 M Q , including 0.20 M Q of 
Fe and 0.02 M of Ni in an inner layer 
as well as ~0.44 M of Fe plus Ni un- 
shocked by the reverse shock and conse- 
quently moving in free expansion. The en- 
ergy of this supernova is estimated from 
the ejecta to be E ps (4-5) xlO 50 ergs 
(Hughes 2000). 

8.2. Li Abundance 

As shown in Figure [6], the Li line at 
6708 A is pronounced in Tycho G. The 
abundance of this element, A (Li) = 2.50 ± 
0.09, raises the question of how much Li 
could have survived in this star that has al- 
ready evolved off the main sequence. Sim- 
ilar abundances have been found in some 
turn-off stars (with 6300 < T cS < 6800 K) 
in the sample of do Nascimento et al. 
(2000, 2003), and in a few subgiant stars 
with stellar parameters similar to those of 
Tycho G. In the sample of Randich et al. 
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Fig. 12. — Li abundance of Tycho G in comparison with that of subgiant stars from Randich 
et al. (1999) having similar stellar parameters, metallicity, and rotational velocity. 



(1999), there are also a few subgiants with 
similar stellar parameters and Li abun- 
dance as Tycho G; see Figure [£21 Al- 
though the Li abundance of Tycho G is 
somewhat high for its relatively low effec- 
tive temperature and surface gravity, the 
upper-right panel of Figure d2] shows two 
stars (with T eff « 5500 K and 6000 K, 
and log (g /cm s2) « 3.3 dex) that are even 
more anomalous than Tycho G. 

Canto Martins et al. (2006) found a Li- 
rich subgiant star (S1242) in the open clus- 
ter M67, with T cS = 5800 K, log(#/cm s2) 
= 3.9 dex, [Fe/H] = -0.05, and A(Li) = 
2.7, all similar to the properties of Tycho 
G. This star is in a wide eccentric binary 



system, and these authors proposed that 
the Li has been preserved due to tidal ef- 
fects. There exist other subgiant stars of 
similar type showing high Li abundances 
and belonging to binary systems. Again, 
this could be linked to being in a close bi- 
nary system where tidal interaction might 
inhibit Li destruction. 

An alternative is that Li was created as 
a result of energetic processes. Martin et 
al. (1992, 1994a,b) suggest that in soft 
X-ray transients (SXTs) one is observing 
freshly synthesized Li, and that the Li 
could come from aa and/or other spalla- 
tion reactions during the outbursts of these 
binary systems that contain black holes or 
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neutron stars. Spallation reactions are pro- 
duced by protons and a particles with en- 
ergies above a few MeV when hitting C, N, 
and O nuclei, and also by a + a collisions 
at similar energies. 

In addition, in the case where the com- 
pact object is a white dwarf, Martin et 
al. (1995) have shown that their com- 
panion stars do not show high Li abun- 
dances. They claimed that this is con- 
sistent with the Li production scenario in 
SXTs, because the weaker gravitational 
potential well of the white dwarfs does not 
allow particles to be accelerated to the re- 
quired high energies (Martin, Spruit, & 
van Paradijs 1994b). However, other au- 
thors have shown that Li production is pos- 
sible during nova explosions (Starrfield et 
al. 1978; Boffin, Paulus, & Arnould 1993). 

The energy arguments could also work 
for the enhancement of Li in the compan- 
ion of a carbon-oxygen white dwarf. The 
accretion lasts for a sufficiently long time 
to build enough Li nuclei in the convective 
envelope of the secondary star. At the on- 
set of the explosion, the supernova is able 
to unbind part of the material from the 
convective envelope. If part of the convec- 
tive envelope were to survive, it would dis- 
play a high Li abundance. However, this 
scenario of Li production during outbursts 
in SXTs also predicts Li isotopic ratios 
in the range 0.1 < iV( 6 Li)/iV( 7 Li) < 10; 
see Casares et al. (2007), and references 
therein. These authors conclude that the 
preservation scenario by tidal effects is the 
most likely mechanism to explain the high 
Li abundances of SXTs. 

There are two other alternatives for en- 
richment of Li in the companions of SNe la. 
The first is that the high-energy particles 



accelerated by the shock wave in the out- 
ermost layers of the exploding white dwarf 
induce spallation reactions in the surface 
of the companion. A fraction of the ir- 
radiated material could mix with the lay- 
ers now making the surface of the compan- 
ion, and Li enrichment would show. The 
second alternative is that the surface of 
the companion, after the explosion, can be 
bombarded by high-energy particles accel- 
erated in the SN remnant, producing spal- 
lation reactions. Steady bombardment by 
those high-energy particles has the restric- 
tion that while it enriches Li, the energy 
deposited would increase the star's lumi- 
nosity. However, if Li enrichment is limited 
to a thin layer extending not much below 
the photosphere, the increase in luminosity 
would be negligible. 

9. Overall Discussion of Tycho G 

Ruiz-Lapuente et al. (2004) suggested 
Tycho G as the mass-donating compan- 
ion of SN 1572 and gave stellar parameters 
that classified it G0-G2 IV star. In 
the present study, higher-resolution spec- 
tra have enabled us to pin down the stel- 
lar parameters: T e fj = 5900 ± 100 K, 
log(#/cm s2) = 3.85 ±0.30 dex, and [Fe/H] 
= —0.05 ± 0.09. Thus, we have been able 
to confirm that Tycho G has a metallic- 
ity close to solar and its luminosity class is 
that of a subgiant near the main sequence. 
Recently, Badenes et al. (2008) have esti- 
mated the global metallicity of the Tycho 
SN la remnant at log(Z SN /Z ) = 0.45±g;|J 
when adopting the solar abundances of 
Grevesse & Sauval (1998). We can esti- 
mate the global metallicity of the star by 
adding all the values of [X/H] provided in 
Tabled] The metallicity of Tycho G would 
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be log(Z TychoG /Z ) a* 0.23, which is con- 
sistent with the value given by Badenes et 
al. (2008) within their large error bars. 

Moreover, in support of the hypothesis 
that Tycho G is physically within the rem- 
nant of SN 1572 and is possibly the surviv- 
ing companion of the supernova, we have 
found a high value of [Ni/Fe], 0.16 ± 0.04, 
about 3<T above the average value in Galac- 
tic disk stars ([Ni/Fe] « -0.05 ± 0.03). In 
addition, we report an anomalously high Li 
abundance for a subgiant, A(Li) = 2.50 ± 
0.09, although similar to a few other known 
subgiant stars having similar stellar pa- 
rameters and metallicity. 

The star is moving at high speed both 
in radial velocity and proper motion. If 
we consider the possibility that Tycho G 
is actually beyond the supernova remnant, 
at 6 kpc or more, where f r ,LSR = —80 
km s~ x would follow the low-latitude rota- 
tion pattern of the Galaxy, we run into two 
other physical incompatibilities: (a) at 6 
kpc in the direction of SN 1572, the re- 
gion has metallicities [Fe/H] < —0.3, much 
lower than observed for Tycho G; and (b) 
the tangential velocity implied by the ob- 
served proper motion and distance would 
be unreasonably high (vt ~ 187 km s _1 at 
d = 6 kpc). On the other hand, accord- 
ing to its space Galactic velocity compo- 
nents, Tycho G could be just a thick-disk 
star passing near the SN la remnant, but 
its high Ni abundance argues against this 
possibility. 

We find that Tycho G has features in 
common with the companions of black 
holes and/or neutron stars that presum- 
ably originated in Type II SNe. These 
companions, which remain bound to the 
black hole or neutron star created by the 



explosion, are main-sequence stars and 
subgiants with a high Li abundance. The 
companions of core-collapse SNe are pol- 
luted with a elements, but Tycho G has 
an overabundance of Ni (and Co at lower 
significance level) with respect to Fe. The 
overabundance is well above the scatter for 
this ratio in Galactic disk stars, suggest- 
ing that Tycho G could have captured the 
low-velocity tail of the SN 1572 ejecta. We 
use SN la yields to model the possible con- 
tamination of the star from the SN ejecta 
and we find reasonable agreement with the 
observed abundances of Tycho G. 

Thus, although the current evidence is 
not fully compelling, it remains consistent 
with Tycho G being the surviving compan- 
ion star of the white dwarf that exploded 
to form SN 1572, as originally proposed by 
Ruiz-Lapuente et al. (2004). 
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Appendices 

A. Low- Resolution Spectra 

We attempted to determine the spectral type of several stars in the Tycho SN 1572 
field. In Figure [13] we display the best fits to the LRIS spectra of Tycho D, E, F, and G 
(see Ruiz-Lapuente et al. 2004 for identifications). The spectra of Tycho G and F were 
dereddened with E (B — V) = 0.76 mag, whereas the spectrum of Tycho D was dereddened 
with E(B — V) = 0.66 mag since this color excess provides a better fit to the LRIS spectrum. 
Based on the closer position of Tycho D and E to the center of the Tycho SN 1572 field, 
we also applied a color excess of E(B — V) = 0.66 mag to Tycho E. However, we know 
from high-resolution spectra that Tycho E is in fact a double-lined spectroscopic binary, 
so this comparison using a low-resolution composite spectrum is meaningless. The spectral 
classifications of most stars differ from those published by Ihara et al. (2007), but are 
consistent with the results of Ruiz-Lapuente et al. (2004). 




4000 5000 6000 7000 

A, (ANGSTROMS) 

Fig. 13. — LRIS spectra of four targets of the Tycho SN 1572 field compared with template 
low-resolution spectra from the Jacoby et al. (1984) library. Note that telluric lines (B band 
near 6900 A; weak lines around 7200 A) were not removed from the templates. The spectra 
of Tycho G and F were dereddened with E(B — V) — 0.76 mag, whereas E(B — V) — 0.66 
mag was used for the spectra of Tycho D and E. The best fit to each spectrum is also 
shown. Note that when Tycho E is observed at high resolution it seems to be a double-lined 
spectroscopic binary, and therefore the fit shown is not reliable. 
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B. Line List with Equivalent Widths for Tycho G 

Here we provide the complete list of line EWs measured from the high-resolution Keck 
spectra of Tycho G. 
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Table 4 

Equivalent Widths of Lines in Tycho G 



Species 


A 


X 




Hj VV \ 


Species 


A 


X 


i°g (gj) 


J±/W \ 










Cm Al 






lev ) 






Li I 


6707.81 


0.00 


0.175 


38.7 ±3.9 


Mnl 


4502.220 


2.920 


-0.490 


60.3 ± 12.3 












Mnl 


5394.670 


0.000 


-3.500 


51.9 ±4.7 


Cl 


5052.160 


7.680 


-1.420 


38.7 ± 7.0 


Mnl 


5399.470 


3.850 


-0.097 


32.0 ±4.7 


Cl 


6587.620 


8.540 


-1.000 


16.5 ±3.8 


Mnl 


5420.360 


2.140 


-1.460 


38.5 ±4.6 


Cl 


7113.170 


8.650 


-0.770 


23.6 ±3.3 


Mnl 


5432.540 


0.000 


-3.620 


36.0 ±4.5 


Cl 


7115.170 


8.640 


-0.930 


32.2 ±3.3 


Mnl 


6021.800 


3.070 


0.030 


91.0 ±4.4 


Cl 


7116.960 


8.650 


-0.910 


25.3 ±3.3 






















Fe I 


5044.220 


2.850 


-2.040 


75.9 ± 7.0 


Ol 


7771.960 


9.110 


2.831 


94.7 ±2.8 


Fe I 


5247.060 


0.090 


-4.930 


58.0 ±5.5 


Ol 


7774.180 


9.110 


2.061 


82.3 ±2.8 


Fe I 


5322.050 


2.280 


-2.900 


47.6 ± 5.1 


Ol 


7775.400 


9.110 


1.256 


68.0 ±2.8 


Fe I 


5806.730 


4.610 


-0.890 


50.3 ±4.4 












Fe I 


5852.220 


4.550 


-1.190 


42.9 ±4.4 


Nal 


5688.220 


2.104 


-0.625 


117.4 ±4.3 


Fe I 


5855.080 


4.610 


-1.530 


12.4 ±4.4 


Nal 


6154.230 


2.102 


-1.607 


34.7 ±3.8 


Fcl 


5856.090 


4.290 


-1.560 


23.5 ±4.4 


Nal 


6160.750 


2.104 


-1.316 


40.1 ±3.8 


Fe I 


6027.060 


4.080 


-1.180 


56.4 ±4.4 












Fe I 


6056.010 


4.730 


-0.500 


52.7 ±4.3 


Mgl 


4730.040 


4.346 


-2.390 


52.3 ±9.4 


Fe I 


6079.010 


4.650 


-1.010 


30.6 ±4.2 


Mgl 


5711.090 


4.346 


-1.706 


109.1 ±4.3 


Fe I 


6151.620 


2.180 


-3.300 


46.0 ±3.8 


Mgl 


6318.720 


5.108 


-1.996 


43.4 ±3.6 


Fel 


6157.730 


4.070 


-1.240 


52.7 ±3.8 


Mgl 


6319.240 


5.108 


-2.179 


29.3 ±3.6 


Fel 


6165.360 


4.140 


-1.500 


45.8 ± 3.8 












Fel 


6180.210 


2.730 


-2.640 


46.1 ±3.7 


All 


6696.030 


3.143 


-1.570 


37.8 ±4.0 


Fe I 


6188.000 


3.940 


-1.630 


43.5 ± 3.7 


All 


6698.670 


3.143 


-1.879 


26.1 ±4.0 


Fe I 


6200.320 


2.610 


-2.400 


73.0 ±3.6 












Fe I 


6226.740 


3.880 


-2.070 


33.4 ±3.6 


Si I 


5665.560 


4.920 


-1.980 


47.0 ±4.3 


Fe I 


6229.240 


2.840 


-2.890 


27.1 ± 3.6 


Si I 


5690.430 


4.930 


-1.790 


52.3 ±4.3 


Fel 


6240.650 


2.220 


-3.290 


42.3 ± 3.6 


Si I 


5701.100 


4.930 


-2.020 


32.2 ±4.3 


Fel 


6265.140 


2.180 


-2.560 


71.6 ± 3.6 


Si I 


5772.140 


5.080 


-1.620 


53.0 ±4.4 


Fel 


6270.230 


2.860 


-2.580 


41.4 ±3.6 


Si I 


5793.090 


4.930 


-1.910 


44.9 ±4.4 


Fe I 


6392.540 


2.280 


-3.930 


11.8 ± 3.7 


Si I 


5948.550 


5.080 


-1.110 


82.6 ±4.5 


Fe I 


6498.940 


0.960 


-4.630 


37.1 ± 3.7 


Si I 


6125.020 


5.610 


-1.520 


38.9 ±3.9 


Fe I 


6608.030 


2.280 


-3.960 


17.9 ± 3.9 


Si I 


6142.490 


5.620 


-1.480 


36.9 ±3.9 


Fe I 


6627.550 


4.550 


-1.480 


24.7 ±3.9 


Si I 


6145.020 


5.610 


-1.400 


35.2 ±3.9 


Fe I 


6703.570 


2.760 


-3.020 


29.6 ±4.0 


Si I 


6155.150 


5.620 


-0.750 


73.2 ±3.8 


Fel 


6710.320 


1.480 


-4.820 


15.4 ±3.9 


Si I 


6244.480 


5.610 


-1.360 


33.8 ±3.6 


Fel 


6725.360 


4.100 


-2.200 


17.7 ±3.9 


Si I 


6721.860 


5.860 


-1.090 


40.0 ±3.9 


Fe I 


6726.670 


4.610 


-1.050 


33.0 ±3.9 












Fe I 


6733.160 


4.640 


-1.430 


25.0 ±3.9 


Si II 


6371.360 


8.120 


-0.050 


42.7 ±3.7 


Fe I 


6750.160 


2.420 


-2.610 


60.2 ± 3.9 












Fe I 


6752.710 


4.640 


-1.230 


33.1 ± 3.9 


si 


6046.020 


7.870 


-0.510 


18.1 ±4.3 












si 


6052.670 


7.870 


-0.630 


12.5 ±4.3 


Fell 


5234.630 


3.220 


-2.230 


98.1 ± 5.5 


si 


6743.440 


7.870 


-1.270 


8.4 ±3.9 


Fe II 


5425.260 


3.200 


-3.160 


66.8 ±4.6 


si 


6757.007 


7.870 


-0.810 


17.8 ±3.9 


Fe II 


5991.380 


3.150 


-3.530 


37.1 ±4.5 












Fell 


6084.110 


3.200 


-3.780 


29.5 ±4.1 


Cal 


5581.970 


2.520 


-0.650 


93.1 ±4.2 


Fell 


6149.250 


3.890 


-2.720 


49.7 ±3.8 


Cal 


5590.120 


2.520 


-0.710 


79.2 ±4.3 


Fell 


6247.560 


3.890 


-2.350 


66.4 ±3.6 


Cal 


5867.570 


2.930 


-1.570 


15.5 ±4.4 


Fell 


6369.460 


2.890 


-4.130 


29.9 ± 3.7 
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Table 4 — Continued 



Species 


A 


X 


log(s/) 


EW A 


Species 


A 


X 


log(ff/) 


EW A 




(A) 


(eV) 




(mA) 




(A) 


(eV) 




(mA) 


Ca I 


6161.290 


2.520 


-1.220 


58.5 ±3.8 


Fell 


6432.690 


2.890 


-3.560 


61.1 ± 3.7 


Cal 


6166.440 


2.520 


-1.120 


65.4 ±3.8 


Fc II 


7479.700 


3.890 


-3.590 


16.7 ±3.0 


Cal 


6169.050 


2.520 


-0.730 


91.5 ±3.7 


Fell 


7711.730 


3.900 


-2.550 


57.9 ± 2.8 


Cal 


6169.560 


2.520 


-0.440 


98.4 ±3.7 












Cal 


6449.820 


2.520 


-0.630 


99.4 ±3.7 


Co I 


4792.860 


3.250 


-0.150 


37.7 ±8.6 


Cal 


6455.600 


2.520 


-1.370 


46.3 ±3.7 


Co I 


5301.040 


1.710 


-1.930 


18.4 ±5.2 


Cal 


6572.800 


0.000 


-4.280 


25.5 ±3.8 


Col 


5483.360 


1.710 


-1.220 


47.9 ±4.3 












Col 


6093.150 


1.740 


-2.340 


13.9 ±4.1 


Sc II 


5239.820 


1.450 


-0.760 


69.2 ±5.5 












Sc II 


5318.360 


1.360 


-1.700 


19.7 ± 5.1 


Ni I 


5082.350 


3.658 


-0.590 


85.4 ±6.9 


Sc II 


5526.820 


1.770 


0.150 


83.8 ±4.2 


Ni I 


5088.540 


3.850 


-1.040 


28.2 ± 6.9 


Sc II 


6245.620 


1.510 


-1.040 


44.5 ±3.6 


Ni I 


5094.420 


3.833 


-1.070 


20.9 ± 6.9 


Sc II 


6300.690 


1.510 


-1.960 


9.9 ±3.6 


Ni I 


5578.720 


1.680 


-2.650 


64.3 ±4.2 


Sc II 


6320.840 


1.500 


-1.840 


12.6 ±3.6 


Ni I 


5587.860 


1.930 


-2.380 


43.8 ±4.3 


Sc II 


6604.600 


1.360 


-1.160 


47.8 ±3.9 


Ni I 


5682.200 


4.100 


-0.390 


59.1 ±4.3 












Ni I 


5694.990 


4.090 


-0.600 


54.8 ±4.3 


Til 


5024.850 


0.818 


-0.560 


52.9 ±7.0 


Ni I 


5805.220 


4.170 


-0.580 


41.9 ±4.4 


Til 


5219.710 


0.021 


-2.240 


20.9 ±5.7 


Ni I 


5847.000 


1.680 


-3.410 


26.4 ±4.4 


Til 


5490.150 


1.460 


-0.980 


23.9 ±4.2 


Ni I 


6086.280 


4.260 


-0.440 


46.3 ±4.1 


Til 


5866.460 


1.070 


-0.840 


36.9 ±4.4 


Ni I 


6111.070 


4.090 


-0.800 


47.9 ±4.0 


Til 


6258.110 


1.440 


-0.440 


47.6 ±3.6 


Ni I 


6128.980 


1.680 


-3.370 


22.2 ± 3.9 


Til 


6261.110 


1.430 


-0.490 


38.4 ±3.6 


Ni I 


6130.140 


4.260 


-0.950 


32.8 ± 3.9 


Til 


6312.240 


1.460 


-1.580 


7.6 ±3.6 


Ni I 


6175.370 


4.089 


-0.550 


56.0 ±3.7 












Ni I 


6176.820 


4.088 


-0.260 


74.0 ±3.7 


vl 


5668.370 


1.080 


-1.000 


5.5 ±4.3 


Ni I 


6177.250 


1.826 


-3.510 


9.4 ±3.7 


vl 


5670.850 


1.080 


-0.460 


9.8 ±4.3 


Ni I 


6204.610 


4.088 


-1.110 


28.7 ±3.6 


vl 


5727.050 


1.080 


-0.000 


39.0 ±4.3 


Ni I 


6378.260 


4.154 


-0.830 


46.4 ±3.7 


vl 


5727 660 


1.050 


-0.890 


15.7 ± 4.3 


Ni I 


6643 640 


1.676 


-2.030 


87 9 + 3 9 


Vl 


5737.070 


1.060 


-0.770 


7.3 ±4.4 


Ni I 


6772.320 


3.658 


-0.970 


51.3 ± 3.8 


Vl 


6111.650 


1.043 


-0.710 


13.1 ±4.0 


Ni I 


7748.890 


3.700 


-0.380 


87.7 ±2.8 


Vl 


6216.350 


0.280 


-0.900 


25.4 ±3.6 


Ni I 


7797.590 


3.900 


-0.350 


73.6 ± 2.8 


Vl 


6251.830 


0.287 


-1.340 


22.9 ±3.6 






















i T 
Zn 1 


4722.160 


4.030 


-0.370 


66.9 ± 9.4 


CrI 


5304.180 


3.460 


-0.680 


14.3 ±5.2 


Zn I 


4810.537 


4.080 


-0.130 


82.3 ± 8.4 


CrI 


5312.860 


3.450 


-0.580 


13.4 ±5.1 


Zn I 


6362.350 


5.790 


0.140 


28.7 ±3.7 


CrI 


5318.770 


3.440 


-0.710 


10.1 ±5.1 












CrI 


5480.510 


3.500 


-0.830 


10.2 ±4.3 


Sri 


4607.340 


0.000 


1.906 


41.4 ±10.8 


CrI 


5574.390 


4.450 


-0.480 


2.7 ±4.2 












CrI 


5783.070 


3.320 


-0.400 


37.4 ±4.4 


Yll 


4900.120 


1.030 


-0.090 


56.5 ± 7.4 


CrI 


5783.870 


3.320 


-0.150 


31.5 ±4.4 


Yll 


5087.430 


1.084 


-0.160 


50.0 ± 6.9 


CrI 


5787.920 


3.320 


-0.110 


42.0 ±4.4 


Yll 


5200.420 


0.992 


-0.570 


51.9 ± 5.8 


CrI 


6330.100 


0.940 


-2.900 


34.3 ±3.6 


Yll 


5402.780 


1.839 


-0.440 


15.2 ±4.7 


Cr II 


5305.870 


3.830 


-1.970 


39.9 ±5.2 


Ball 


5853.690 


0.604 


-0.910 


78.7 ±4.4 












Ball 


6141.730 


0.704 


-0.030 


134.9 ± 3.9 


Mnl 


4265.920 


2.940 


-0.440 


55.8 ± 13.5 


Ba II 


6496.910 


0.604 


-0.410 


112.6 ± 3.7 


Mn I 


4470.130 


2.940 


-0.550 


56.3 ± 12.4 
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